GMOS-IFU spectroscopy of 167-317 (LV2) proplyd in Orion 1 
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ABSTRACT 

We present high spatial resolution spectroscopic observations of the proplyd 
167-317 (LV2) near the Trapezium cluster in the Orion nebula, obtained during 
the System Verification run of the Gemini Multi Object Spectrograph (GMOS) 
Integral Field Unit (IFU) at the Gemini South Observatory. We have detected 
38 forbidden and permitted emission lines associated with the proplyd and its 
redshifted jet. We have been able to detect three velocity components in the 
profiles of some of these lines: a peak with a 28-33 km s _1 systemic velocity that is 
associated with the photoevaporated proplyd flow, a highly redshifted component 
associated with a previously reported jet (which has receding velocities of about 
80-120 km s" 1 with respect to the systemic velocity and is spatially distributed 
to the southeast of the proplyd) and a less obvious, approaching structure, which 



1 Based on observations obtained at the Gemini Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the 
Gemini partnership: the National Science Foundation (United States), the Particle Physics and Astronomy 
Research Council (United Kingdom), the National Research Council (Canada), CONICYT (Chile), the 
Australian Research Council (Australia), CNPq (Brazil) and CONICET (Argentina). 

2 Current Address: ICN-UNAM, Ap. Postal 70-543, CU, D.F. 04510, Mexico 
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may possibly be associated with a faint counter-jet with systemic velocity of 
(—75 ± 15) km s -1 . We find evidences that the redshifted jet has a variable 
velocity, with slow fluctuations as a function of the distance from the proplyd. 
We present several background subtracted, spatially distributed emission line 
maps and we use this information to obtain the dynamical characteristics over 
the observed field. Using a simple model and with the extinction corrected Ha 
fluxes, we estimate the mass loss rate for both the proplyd photoevaporated flow 
and the redshifted microjet, obtaining M prop iyd = (6.2 ± 0.6) x 1CT 7 M year -1 
and M jet = (2.0 ± 0.7) x 10~ 8 M year -1 , respectively. 

Subject headings: ISM: clouds — ISM: individual (Orion Nebula, LV2, 167-317) 

1. Introduction 

The Orion Nebula (M42) is the most active site of star formation of the Orion Molecular 
Cloud and contains many young stars, among them the Trapezium cluster. Its high mass 
young stellar members generate an intense ultraviolet radiation field that photodissociates 
and photoionizes the nearby material. These mechanisms promote the appearance of some 
complex structures as the so-called proplyds (O'Dell & Wen 1994). Proplyds are low mass 
YSOs that are being exposed to an intense ultraviolet radiation field which renders them 
visible. In Orion, there are approximately 160 proplyds (Bally et al. 2000) which are being 
photoionized mainly by 9 1 Ori C, an 06 spectral type star. There are several objects that 
have been identified as being proplyds not only near the Trapezium cluster (Bally et al. 2000) 
but also in other star forming regions (like NGC 3372, Smith et al. 2003). Most of them 
share the same features: a bow-shaped head that faces the ionization source, a tail, that is 
primordially directed away from the source, a young star, that may (or not) be visible, and 
a disk, sometimes seen in silhouette against the HII region (e.g., O'Dell et al. 1993; O'Dell 
1998; Bally et al. 2000; Smith et al. 2005). Early studies of these objects, however, were only 
able to determine their apparent ubiquity close to the 9 l Ori C star, their high ionization 
level and, in later studies, the presence of high velocity structures (outflows) associated with 
these condensations (e.g., Laques & Vidal 1979; Meaburn 1988; Meaburn et al. 1993). 

The proplyds are presently explained by a set of models which include a photoevapo- 
rated wind, an ionization front and a photoionized wind (Johnstone et al. 1998; Storzer & 
Hollenbach 1998; Henney & Arthur 1998; Bally et al. 1998; O'Dell 1998; Richling & Yorke 
2000). Johnstone et al. (1998) proposed models in which far ultraviolet radiation (FUV) and 
extreme ultraviolet radiation (EUV) are responsible for the mass loss rates of the proplyds 
depending on the distance of the proplyd to the ionization source. For the proplyds situated 
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at intermediate distances from the ionizing star, FUV photons penetrate the ionization front 
and photodissociate and photoevaporate material from the accretion disk surface, generating 
a supersonic neutral wind that passes through a shock front before reaching the ionization 
front as a neutral, subsonic wind. At the ionization front, EUV photons ionize the wind, and 
the material is then reacelerated to supersonic velocities. The interaction of this supersonic 
ionized wind with the star wind generates bow shocks seen in Ha and in [O III] A5007 in some 
proplyds (Bally et al. 1998). For proplyds closer to the ionizing star, the ionization front 
reaches the disk surface, and the resulting wind is initially subsonic, becoming supersonic 
further out. 

In general, the models cited above are able to explain the main observational features 
of the Orion proplyds. The region close to the ionization front is responsible for most of 
the emission of these objects. Storzer & Hollenbach (1998) have shown that heating and 
dissociation of H 2 can explain the observed [O I]A6300 emission at the disk surface (Bally et 
al. 1998). Bally et al. (2000) proposed that neutral gas compressed between the shock front 
and the ionization front can explain the [O III] A5007 emission seen near the ionization front 
in some objects. Henney & Arthur (1998) were able to reproduce the observed Ha intensity 
profile of the proplyds in Orion in terms of accelerating photoevaporated flows. Numerical 
2D HD simulations with a treatment of the radiative transfer (Richling & Yorke 2000) also 
reproduce the proplyd morphology and emission, although in this work the treatment of the 
diffuse radiation is somewhat simplified. 

There are several open issues about these objects. One of them is related with the 
calculation of the mass loss rate of the proplyds, that is strongly model dependent and 
which may pose severe constraints on the age of the Orion proplyds (as well as for 9 1 Ori 
C) . As the flux of FUV photons dissociates the molecules in the disk of the proplyds, the 
sub- or trans-sonic wind is accelerated at the proplyd's ionization front (IF) (see Henney et 
al. 2002 for a discussion), leading to derived mass- loss rates of 8 x 10~ 7 M Q yr _1 which 
imply a short lifetime for these systems (Churchwell et al. 1987). It is hard to reconcile 
such short lifetimes with the fact that the region of the Trapezium cluster is populated by 
several proplyds. Improvements in both models and observational techniques have slightly 
reduced the calculated mass-loss rates for these systems, and Henney et al. (2002) find 
M = 8.2 x 1O~ 7 M yr" 1 for the 167-317 1 (LV2) proplyd, implying an age of less than 10 5 
yr for 6 1 Ori C, and M < 1.5 x 1O- 6 M yr" 1 for the 170-377, 177-341, 182-413, 244-440 
proplyds (e.g., Henney & O'Dell 1999). 



1 Thc proplyds mentioned in the present paper will be denoted following the O'Dell & Wen (1994) notation, 
that is based on the coordinate of the object. 
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It is presently known that many proplyds show, besides protostellar features such as 
accretion disks and a young low mass protostar, the presence of jets. Bally et al. (2000), in 
a survey carried out using the WFPC2 camera on the HST, find 23 objects which appear to 
have collimated outflows seen as one-sided jets or bipolar chains of bow shocks. Meaburn 
et al. (2002) also found kinematical traces of jets associated with LV5 (158-323) and GMR 
15 (161-307). The proplyd 167-317 shows evidence of the existence of a collimated outflow, 
detected spectroscopically (Meaburn 1988; Meaburn et al. 1993; Henney et al. 2002). Until 
now, this outflow was detected as a one-sided jet, with a P.A. of 120° and a propagation 
velocity of about 100 km s _1 . Recently, Smith et al. (2005) and Bally et al. (2005), using 
the Wide Field Camera of the Advanced Camera for Surveys (ACS/WFC) on HST, found 
jets associated with silhouette disks in the outer regions of Orion. 

In this work, we present the first Integral Field Unity (IFU) observation of a proplyd. 
The observed object is 167-317, one of the brightest proplyds of the Trapezium. We discuss 
one of the very first results from the Gemini Multi-Object Spectrograph, in its IFU mode 
(hereafter, GMOS-IFU), at the Gemini South Telescope. In §2, we present the observations 
and the steps to obtain a calibrated, clean datacube. In §3 we present the spectral analysis 
for the observed region of the 167-317 proplyd, as well as the observed lines and intensity 
maps. We use the observed line profiles to identify the outflows in the system. In §4, we 
present the discussion and the conclusions. 



2. The observations and data reduction 

The data were taken during the System Verification run of the GMOS-IFU at the Gemini 
South Telescope (GST), under the GS-2003B-SV-212 program, on 2004 February 26 t?i and 
27 th . The science field of view (FOV) is 3"5 x 5" with an array of 1000 lenslets of 0"2. 
The sky is sampled with 500 lenses which are located 1' from the science FOV. We used 
the R831_G5322 grating in single slit mode, giving a sampling of 0.34A per pixel 15 km 
s _1 per pixel at Ha) and a spectral coverage between ~ 5515 A to ~ 7630A. The spectral 
resolution (i.e., the instrumental profile) is 47 km s" 1 < FWHM < 63 km s" 1 (the FWHM 
decreasing with increasing wavelength). A field centered on the 167-317 proplyd has been 
observed with an exposure time of 300s. A 60 s exposure of the standard star the Hiltner 
600 has also been taken in order to derive the sensitivity function. 

The data were reduced using the standard Gemini IRAF vl.6 2 routines. An average 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
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bias image has been prepared using the GBIAS task. The extraction has been worked out 
using the GCAL flat and the response using the twilight flat. Spectra extraction has been 
performed and wavelength calibration done using an arc taken during the run. The 167-317 
field and Hiltner 600 star have been processed in the same way. The sensitivity function 
derived from the standard star has been applied to the observed field. Finally, data cubes 
have been built using the GFCUBE routine. We maintain the IFU original resolution of / /2 
px _1 although an interpolation was performed in order to turn the IFU hexagonal lens shape 
into squared pixels. 

Because some of the emission lines are very intense, it was not possible to eliminate 
cosmic rays using the standard GSCRREJ routine of the IRAF Gemini vl.6 package. In order 
to remove cosmic rays we have then programmed an IDL routine. This routine works directly 
on the data cube by taking out impacts that are above a certain convenient level calculated 
based on the mean intensity of nearby pixels. The cosmic ray is then eliminated doing a 
linear interpolation in the wavelength direction. With this procedure, it is possible to keep 
very intense lines and eliminate lower level cosmic ray impacts that could be misinterpreted 
as low emission lines. 

Further manipulation of the cube, such as the production of channel maps, emission 
line maps, dispersion and velocity maps, have been done using Starlink 3 , IDL routines and 
Fortran programs developed specially for this purpose. 

The subtraction of the background emission from the proplyd spectrum is a very chal- 
lenging task (Henney & O'Dell 1999). This problem arises since the background nebular 
emission is strongly inhomogeneous on an arc-second scale. In this work, we have used 
standard x 2 - m ethod in order to fit an intensity versus position plane in a semi-rectangular 
field located ~ 1" — l'/5 away from the peak proplyd emission in each of the individual ve- 
locity channel maps. We assure that we only take samples of the background avoiding the 
region of the proplyd. These planar fits to the nebular emission are then subtracted from 
the corresponding velocity channel maps. In Figure 1, we show the observed field, the object 
orientation in the plane of the sky, and also the region that defines our background (region 
4, labeled as R4 in Figure 1). Region 4 (see Figure 1) contains 76 spectra which are used 
to define the background. This Figure also shows an internal box (labeled as Rl, or region 
1), for which we have defined an arbitrary xy-coordinate system limited by 1 < x < 11 (in 
the East direction) and 1 < y < 17 (in the North direction; see Figure 1). In the following 
sections, we show maps of spatially distributed physical variables that will be constrained 



Science Foundation. 

3 See documentation at this website: http://star-www.rl.ac.uk/ 
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to the limits defined by this box. The spatial positions inside this domain will be defined 
using the xy-coordinate system described above. Regions 2 and 3 (labelled as R2 and R3 in 
Figure 1) define, respectively, a region near the center of the 167-317 proplyd and a region 
which has a high- velocity, redshifted feature (see below). 

3. Observational results: spectral line identification and high velocity features 

3.1. The emission lines intensities and ratios 

The lines that we identify in the spectra have already been reported in previous papers 
of the Orion Nebula (e.g,. Baldwin et al. 2000). We select 4 different regions to measure the 
line intensities, namely, regions Rl, R2, R3 and R4 (see Fig. 1 and §2 for the coordinate 
system definition). 

Four spectra are then obtained by co-adding the spectra of the pixels included in each 
region. The spectrum of region 1 is the result of the sum of 187 pixels, for region 2 the 
sum of 6 pixels, region 3 represents the spectrum of a single pixel and region 4 the sum of 
76 pixels. A set of 38 lines has been selected from a manual search for emission features 
in the spectrum integrated over the selected fields. The intensities of the emission lines are 
determined by fitting a flat continuum and integrating over the whole emission feature 4 . 
Then, a mean flux was obtained dividing the integrated flux by the total number of pixels of 
each region. Table 1 gives the list of the observed lines together with the mean flux of each 
region (after subtraction of the emission from the background, region R4) and the ratio to 
Ha (normalized to Ha = 100). The fourth and eigth columns show the background emission 
and ratios to Ha, respectively, for each line. The mean fluxes are given in units of 10~ 15 
erg cm~ 2 s _1 px _1 , and are not corrected for reddening. The absolute error of the mean 
flux is given in parentheses. This error was calculated taking into account the variations of 
the local continuum. The errors can be very large for the weaker lines. Also, the detection 
of spectral lines in region R3 is affected by the low signal to noise ratio. There are some 
lines which show stronger mean fluxes in the background, namely, the lines of Si III A5740, 
Si II A5979, [Ni II] A7378 and [Fe II] A7155 (see Table 1). These lines have low ionization 
potentials and are expected to appear mainly in regions of lower ionization degree. The 
lower line intensities within the proplyd indicate that these lines are mainly emitted by the 
background nebula, and that they are absorbed at least partially by the dust in the proplyd. 



4 The radial velocity range that defines the integration limits in each intensity line determination varies 
from line to line, since the FWHM changes with the wavelength (see §2). 
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In Figure 2, we show the full extracted spectra from the data cube, for regions Rl to R4 
as defined in Figure 1. In this figure, the spectra for regions Rl, R2, R3 and R4 are depicted 
from bottom to top in each diagram. The data shown here are not background subtracted. 
Most of the detected 38 lines can be clearly seen in the spectra. It is also evident that the 
S/N is higher in regions Rl and R4 (which is due to the increase of the pixel number of these 
regions). Here, a third order cubic spline was fitted to the continuum and then subtracted 
from the data. Because of the high order spline polynome subtraction, some minor variations 
are still present near the Ha line. 

Figure 3 depicts (for region Rl), 2D-intensity maps of the lines [N II] A5755, [N II]A6548, 
Ha, [N II] A6583, [He I]A7065 and [Ar III] A7135, superimposed by line profiles. The data 
shown here are background subtracted, as described in §2. The line profiles can be better 
seen in Figure 4, in which we show the mean line profile for each line and for each of the 
regions Rl, R2, R3 and R4, in the same diagram. As could be anticipated, the flux is more 
intense in region R2 (dotted line), where the proplyd is located, except for the [N II]A6548 
line, which shows a stronger flux in region R3. For the Ha and [N II] A6583 lines the emission 
of region R3 is more intense than the emission of region Rl. For lines with good S/N, the 
presence of a redshifted feature starts to become clear. An analisys of the profile components 
will be given in more detail below. 

In order to see how the spectra of the different lines change spatially, in Figure 5 we 
present an intensity versus position plot for the same lines shown in Figures 3 and 4. We 
show how the intensity of these lines changes along a diagonal line, that crosses region Rl 
passing through the pixel (3,4) (region R3) and through the centre of region R2. We note 
that there are 2 intensity maxima for the Ha, [N II] A6583 and [N II]A6548 lines, clearly 
showing in which lines region R3 is clearly visible (also see Figure 3). 

3.2. High velocity features 

We associate region R3 with the redshifted jet of the proplyd, although in Figure 9 we 
will show that the emission of the jet is not restricted only to this pixel. A redshifted jet 
associated with the 167-317 proplyd was previously detected, with a propagation velocity of 
100 km s _1 (e.g., Henney et al. 2002). We have also detected this high velocity feature in 
our data. In order to see this, we examine the behaviour of the moments of the radial velocity 
distribution. These are the flux- weighted mean radial velocity < v >, the flux- weighted rms 
width of the line Av 2 and the skewness Av 3 , that are given, respectively, by: 
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The integrated intensity (see eq. 4 and Figure 3) shows that if we move away from the 
prop lyd position (R2), the emission for several lines (for example, Ha and [Ar III]) rapidly 
drops to low values (close to the background value). We also note that the maps for the 
flux-weighted mean radial velocity 5 (not shown here; see eq. 1) show that < v >= 50 ± 20 
km s _1 . 

The flux- weighted rms width of the line (At; 2 ) 1 / 2 and the skewness (Aw 3 ) 1 / 3 moment of 
the radial velocity distribution (see equations 2 and 3, respectively) are depicted in Figure 
6, for the Ha (left) and [Ar III] (right) lines. Both of them were computed after carrying out 
both background and continuum subtractions. From Figure 6, we see that there is a clear 
enhancement of both of these moments from the proplyd position (R2 region) towards the 
R3 region (i.e., towards the SE direction, at a PA ps 135°). This behaviour of the moments 
leads us to infer that: 1) there is an increase in the line width as we go from the proplyd 
position to the SE direction; and 2) that there are redshifted wings in the line profiles. These 
results indicate the presence of a redshifted outflow in this region. We note that the same 
behaviour is seen in the He I A6678.15, [N II]A6583.46 lines (not shown here). 

In order to confirm the presence of high velocity components in the observed field, as well 
as to see the behaviour of the photoevaporated proplyd flow and the background emission, 
we have computed three-component Gaussian minimum x 2 fits for each position-dependent 
line profile. The profiles of all emission lines have a major peak at v ra d = 40 — 60 km s _1 , 



5 Unless explicitly mentioned, the radial velocities presented here are not corrected for the systemic radial 
velocity, that is of the order of v Q w 26 km s -1 (e.g., Meaburn et al. 1993). The error in the radial velocity 
measurements is of the order of 2.5 km s _1 . 
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the exact radial velocity of this peak changing with spatial position. In some spectral lines, 
there is an evident second peak at redshifted velocities ranging from 100 to 150 km s -1 , 
corresponding to the jet associated with the proplyd. Finally, we have been able to detect 
a blueshifted component, which is fainter than the redshifted components. As an example 
of our three-component Gaussian fit, in Figure 7 we show the data (full line) and the fit 
(crosses), for the Ha and [Ar III] A7135 lines in the (x,y)—(6,7) position. In this figure, 
for each emission line, the top-left panel represents the data and the three Gaussian fits, 
the top-right panel shows the main, low velocity component (and the Gaussian fit for this 
component), the bottom- left panel depicts the data minus the main component together with 
the fits for the blue- and red-shifted components and, finally, the bottom-right panel shows 
the residual, obtained by subtracting the three-component fit from the observed line profile. 
The fits depicted in this figure show the presence of redshifted and blueshifted components. 
The redshifted component is present in several lines in pixels around the SE direction, as 
already mentioned before. On the other hand, blueshifted emission can be found in several 
pixels to the NW of the proplyd peak emission as can be seen in the left bottom panel of 
Figure 7 and this could represent the first spectroscopic determination of the presence of a 
blueshifted counter-jet, with systemic velocities of -60 km s _1 < v rad < —90 km s" 1 . This 
blueshifted component is also suggested by the fitted profiles of Hel A7065 (not shown here). 

Figure 8 depicts the spatially distributed intensity over the Rl region (see Figure 1), of 
the main peak for the [N II]A5754 (top- left), [N II]A6548 (top-middle), Ha (top-right), [N 
II]A6583 (bottom-left), Hel A7035 (bottom-middle) and [Ar IIIJA7135 (bottom-right) lines. 
The values shown in this Figure were obtained from three-component Gaussian fits to the 
observed line profile. All of the lines peak at the proplyd position, at (x,y) = (6,7), inside 
the R2 region, and three of them (namely [N II]A6548, top-middle; Ha, top-right; and [N 
II] A6583, bottom-left), have a secondary, less intense peak at (x,y) = (3,4) (i.e., in the R3 
region; see Figure 1 for the definition of the coordinate system). This secondary peak is 
probably related to the jet, since, as we have seen before, the jet propagates from the R2 
region towards the SE direction. It is also interesting to note that there is a tail of faint 
emission (compared with the maximum in a given map) that extends in the NE direction, 
and that can be seen in both of the [N II] lines that bracket Ha (top-middle and bottom-left 
maps in Figure 8). 

In Figure 9, we show the spatially dependent intensity (Fig. 9a) and the central velocity 
(Fig. 9b) of the redshifted component of the fitted Ha profile. We limit the maps to the 
spatial pixels which have a skewness | (At; 3 ) 1 / 3 1 > 30 km s _1 . The pixels which satisfy this 
criterium have a well defined high velocity, redshifted component. These figures show that 
the high velocity component has an intensity maximum near the center of the proplyd, 
with an extension towards the SE direction (along the jet axis), surrounded by a region of 



-10- 



decreasing fluxes. The high intensity spike seen in Ha, [N II]A6548 and [N II] A6583 lines 
(Figure 3) can be seen here. The jet shows a spread in velocity values, although most of the 
pixels present velocities around 120-140 km s _1 . This figure also shows that the jet velocity 
decreases with increasing distance from the proplyd. This trend can be seen in the Figure 
9c, where we show the radial velocity of the redshifted component as a function of distance 
from the proplyd position, (x,y) = (6,7) in region Rl. To build this figure, we assume that 
the jet is propagating in a PA 135° position angle. We define a box with one axis aligned 
with the jet axis, and the second, perpendicular dimension extending two pixels to each side 
of the jet axis. We then take averages (of the mean velocity of the redshifted component) 
perpendicular to the jet direction and plot the resulting velocity as a function of distance 
from the proplyd (see the bottom frame of Figure 8). There is an indication that the jet 
velocity is slowly diminishing as a function of distance from the source. 

3.3. The spatial distribution of the nitrogen ratio 

In Figure 10a, we show a map for the Rl region of the line ratio: 



[N HJA6548 + [N II]A6583 = A 654 8^6548 + ^6583^6583 (5] 

[N II]A5754 A 5754 (J 5754 ' 1 ' 

where A and a are the height and the dispersion of the fitted Gaussian profile (of the 
main, low velocity component). The [Nil] (6548+6583) /5754 ratio is a classical electron 
temperature diagnostic of low/medium density nebulae (N e < 10 5 cm -3 ). At high densities, 
due to the collisional desexcitations, this ratio cannot be used for the determination of the 
electron temperature, but it turns out to be a quite good electron density diagnostic. This 
can be seen, for example, from the Osterbrock (1989) equation: 



6.91exp(2.5 x 10 4 /T e 
1 + 2.5 x W- 3 (N e /Te /z ) 



[Nil] (6548 + 6583)/5754 = ^ !^% (6) 



Figure 10b shows the lines corresponding to [Nil] (6548+6583) /5754 = 9, 10, 11, 80, 90, 
100 in a T e — N e plane. We can see that for the values obtained in our case (line ratios close 
to 10), the [Nil] (6548+6583) /5754 ratio is strongly dependent on the electron density, and 
that the dependence on the electron temperature is small. 

In Figures 10c and lOd, we show the densities computed from the [NII](6548+6583)/5754 
ratio for T = 10 4 K and T = 1.5 x 10 4 K, respectively. Although presenting a very com- 
plicated pattern, these figures show that the region coinciding with the proplyd center, and 
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extending towards the SE direction (the jet propagation direction), shows the lowest line 
ratio, indicating electron densities (n e ) of the order of 2 x 10 6 cm -3 . This value is at least 
in qualitatively agreement with the no = (3.0 ± 0.5) x 10 6 cm -3 obtained by Henney et al. 
(2002) considering a T = 1.2 x 10 4 K temperature for the ionization front (that is inside the 
R2 region). We can also note that both figures are similar, presenting few differences in the 
density values. 

The extended structure from the R2 region towards the SE direction is surrounded by 
a region in which the nitrogen line ratio reaches values of up to w 30. In this limit, n e < 10 5 
cm~ 3 . It is interesting to note that in the clump associated with the redshifted jet, the 
nitrogen line ratio increases substantially, indicating a decrease in the electron density when 
compared with the same values near the center of the proplyd. We have also obtained the 
[S II] A6716, A6730 emission maps. Unfortunately, we are not able to use these lines as a 
diagnostic because neither the ratio of this doublet is a good electron density indicator for 
this case (since it is constant for densities higher than 10 5 cm~ 3 ) nor do we have a good 
enough signal to noise ratio for these lines. Furthermore, the background subtracted spectra 
for these two lines give us negative fluxes in the proplyd region (see Table 1). The presence 
of these negative features in the line profiles is probably related with the presence of dust, 
as discussed by Henney & O'Dell (1999). 



4. Discussion and conclusions 

We have presented in this work the first Integral Field Unit (IFU) spectroscopic obser- 
vations of a photoevaporating disk immersed in a HII region. In particular, we have taken 
advantage of the System Verification Run of the Gemini South Telescope Multi Object Spec- 
trograph (GMOS) to obtain spectra of the 167-317 proplyd in the Orion nebula. The 167-317 
proplyd, also known as LV2 (from the pioneering work from Laques & Vidal 1979), is one 
of the brightest and best studied proplyds. In a single exposure, we took 400 spectra with 
a spatial resolution of / /2. These spectra have been combined in order to optimize this 
instrumental feature, and, before discussing of our results, we want to make a few comments 
concerning the potential use of such IFU data-cubes for these objects. As previously dis- 
cussed in the literature (Henney et al. 2002), the subtraction of the background emission is a 
very challenging task. This problem arises since the background near these proplyds (which 
are immersed in a highly non- homogeneous, photoionized ambient medium) is very complex 
and changes in an arc-second spatial scale. A careful computation of the background con- 
tribution to the spectra of these objects is one of the most important tasks in order to have 
reliable radial velocity measurements for both the proplyd photoevaporated flow and the 
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high velocity features (jets) present in these systems. Here, we have used the 76 spectra, 
defined in Figure 1 as region R4, in order to obtain a planar fit for the background emission. 
We think that such a definition of the background emission is more precise, since we use 
several regions to calculate a better function to describe the background. One more thing 
to point out is that, with spatially distributed spectra, we are able to confidently separate 
the contribution of the background from the contribution of the object itself. 

The data from the IFU observation of the 167-317 proplyd also allow us to investigate 
the spatially dependent properties of the outflows associated with this object. In particular, 
the 167-317 proplyd is known to have a redshifted, collimated jet that propagates with 
heliocentric velocities of ~ 100 km s _1 towards the SE direction, and with a spatial extension 
of ~ 2" (e.g, Meaburn 1988; Meaburn et al. 1993; Massey & Meaburn 1995; Henney 2000; 
Henney et al. 2002). We find that a prominent, high velocity redshifted component can be 
detected in some emission lines, particularly in Ha, in this SE direction. The redshifted jet 
has a trend in its radial velocity (as a function of distance from the proplyd) with higher 
velocities close to the proplyd and lower velocities at increasing distances. We note that 
Henney et al. (2002) had previously suggested a variation in the jet velocity as a function of 
distance from the source. 

There is a subtle peak (in several emission lines, particularly in [N II]A5754 and Ha) that 
is associated with the SE jet emission, and located at ~ 0'.'67 SE of the center of the proplyd 
[at (x,y) = (3,4); the R3 region in Figure 1]. For this intensity peak, and also using the Ha 
profiles of the neighbouring pixels, we find a mean heliocentric velocity v re( i = 116 ± 10 km 
s -1 . It is interesting to note that previous HST images and spectroscopic analyses (Bally et 
al. 2000; Henney et al. 2002) also reveal a spike in the jet emission at « 0'/4 from the proplyd 
cusp. If we associate these two emission regions as belonging to the same condensation, we 
can estimate a proper motion of ~ 140 km s -1 , which combined with the ~ 116 km s~ 1 
radial velocity give a full jet velocity of xs 180 ± 90 km s _1 (where the error was inferred 
from the spatial resolution of our data). This velocity is of the same order as the values 
inferred for Herbig-Haro jets associated with T Tauri stars. However, the exact nature of 
this spike is unknown. We could relate this spike with a bow shock of the jet, since the line 
profiles for both are similar (Hartigan et al. 1987; Beck et al. 2004), showing a low velocity, 
more intense peak together with a high velocity, less intense peak (similar to the line profiles 
seen in Figures 3, 4 and 7). 

We have found evidence for a blueshifted component in several emission lines around the 
northwest tip of the proplyd position. In particular, the Ha profiles reveal the presence of 
a blue-shifted component with systemic velocity vuue = (—75 ± 15) km s _1 , which indicates 
that its velocity is similar though smaller than the one of the redshifted jet. However, the 
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intensities of the blue-shifted features are much lower than the intensity of the redshifted 
components associated with the SE jet: the counter-jet is at least four times less intense 
than the redshifted jet (in Ha). The presence of such a blue-shifted component very close 
to the LV2 proplyd, together with previously reported evidence that (at several arcseconds 
to the NW of this proplyd; see Massey & Meaburn 1995) there is a blueshifted component 
(in the [O III] A5007 line) strongly suggests the existence of a faint counterjet. 

It is interesting to estimate the mass loss rate for the proplyd and the associated red- 
shifted jet. In order to do this, we consider a simple model that assumes that the proplyd 
flow arises from a hemispherical, constant velocity wind that originates at the proplyd ion- 
ization front, at a radius tq. Using the Ha luminosity Ln a , we can then obtain the particle 
density n at this point and the mass loss rate, from: 



where e THa accounts for the extintion correction, an a = 5.83 x 10~ 14 cm 3 is the effective 
recombination coefficient for Ha, hv Ha is the energy of the Ha transition, /i = 1.3 is the 
mean molecular weight and c/j = 10 km s -1 is a typical value for the sound velocity in an 
HII region. The extintion correction is obtained from the base 10 logarithm of the extintion 
at HP {c H f}) which for the proplyd 167-317 is equal to 0.83 (O'Dell 1998) through the 
relation, cap = Kr Ha (O'Dell et al. 1992), where here we take K = 0.56. Assuming for the 
ionization front radius a value ro = (7.9 ± 0.2) x 10 14 cm (Henney et al. 2002) we obtain 
n = (2.3 ± 0.6) x 10 6 cm -3 and M = (6.2 ± 0.6) x 10~ 7 M year" 1 . Henney et al. (2002) 
have obtained a M = 8.2 x 10~ 7 M year -1 ±10%, and the small differences between both 
results may arise from the different observational techniques employed in both cases. 

For the redshifted jet, we assume that the outflow arises from the blob located in region 
R3 (or pixel 3,4). In this case, the particle density, mass and mass loss rate of the jet are 
given by, 




(7) 



M 



4irrln iim H cii, 



(8) 




1/2 



(9) 



M = —r s B nfj,m H , 



(10) 
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M ~ My-, (11) 

where tb = 6.7 x 10 14 cm is the radius of the blob, Vj = 180 km s -1 and Lj = 4.5 x 10 15 cm 
are the propagation velocity and length of the jet, respectively, derived above. With these 
values, we obtain for the jet a mass loss rate equal to (2.0 ± 0.7) x 10 -8 M year -1 , which 
is similar to the mass loss rate of typical HH jets from T Tauri stars. 

We have been able to construct a ([NII]A6548 + [NII]A6583)/[NII]A5754 line ratio map. 
This ratio indicates densities higher than 10 5 cm -3 for the proplyd emission region, consistent 
with values previously reported in the literature. We find a subtle enhancement of this ratio 
in the region of the redshifted jet. Further observations of electron temperature and density 
diagnostic lines would allow a more accurate determination of the spatial dependence of 
these important parameters. 
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Table 1: Detected emission lines for 167-317 
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The terms inside the parentheses are the absolute errors of the intensity of each line, for each one of the 
regions. 

The symbol ND implies that the line were not detected in that region. 
Absolute error smaller than 6 x 1CT 4 . 
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Fig. 1. — An Ha integrated emission map (shown with a linear grey scale in arbitrary 
units) showing the location of the 167-317 proplyd with respect to the scientific field of view 
(SFOV), and its orientation: North is up and East is left, as indicated by the arrows on the 
top-left corner of the figure. The field has 16 pixels in the East direction and 25 in the North 
direction. The R4 frame located at the left and top parts of the figure represents the region 
choosen to define the background, which covers two columns of pixels on both the North 
and East directions. An internal box (Rl), approximately centered on the object emission 
defines a coordinate system, and the numbers on the bottom axis and on the left axis of this 
internal box are referred to in the text as the x- and ^/-coordinates, respectively. This spatial 
sub-sample of spectra will be used throughout the paper. The 9 1 Ori C star is located at 
7'.'8 from the proplyd, at a P.A. = 236°. Four regions are defined in this figure, namely Rl 
(the internal box), R2 (the proplyd center), R3 (a clump associated with a redshifted jet) 
and R4 (the background). We will refer to these regions throughout the paper. 



-18- 



1 .c 

it 0.8 

~ E 0.6 

u 

» 0.4 

b 

5- 0.2 



0.0 



1 1 


1 1 1 1 1 


1 1 


1 . . . . 






1 1 1 1 



5500 5600 5700 5800 5900 6000 
A (A) 



.< 0.8 - 



0.6 



0.0 




> 0.4 >MVf»M**W** |, ^')*--MW(,Wf' 



0.2 



6000 6100 6200 6300 6400 6500 
A (A) 




6500 6600 6700 6800 6900 7000 7000 7100 7200 7300 7400 7500 

A (A) A (A) 



Fig. 2. — In each diagram, from bottom to top, are depicted, respectively, the spectra for 
regions Rl, R2, R3 and R4. A 3rd order cubic spline was fitted in order to subtract the 
continuum and a constant n has been added to the spectra of regions R2 (n = 0.2), R3 
(n = 0.3) and R4 (n = 0.4) in order to produce an intensity offset between the successive 
spectra. Wavelengths (A) are in A, and the fluxes are in units of 10 -15 erg cur 2 s _1 A -1 . 
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Fig. 3. — A sub-sample of the observed FOV, namely, the Rl region (see Figure 1), showing 
the integrated flux (the grey-scale map) for the following lines: [N II] A5755 (top-left), [N 
II]A6548 (top-center), Ha (top-right), [NII]A6583 (bottom-left), [He I]A7065 (bottom-center) 
and [Ar III] A7135 (bottom-right). The linear bars give the grey-scale in unit of 1 x 10~ 15 
erg cm -2 s _1 px _1 . The spectra in each position was background subtracted. Superimposed 
on the intensity maps are the line profiles for each box. See the text for a discussion. The 
orientation is the same as in Figure 1 (i.e., North is up and East is left). The original pixels 
(see Figure 1) have been resampled in order to smooth the emission features, and each box 
in these panels are slightly larger than the original pixels in Fig. 1. 
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Fig. 4. — An alternative view of the line profiles depicted in Figure 3. Each diagram shows 
the line profile for regions Rl (continuous line), R2 (dotted line), R3 (dashed line) and R4 
(dash-dotted line). The emission line is indicated on the top-left corner of each plot. The 
fluxes are in units of 1 x 10~ 15 erg cm~ 2 s _1 per A and tey are plotted against the radial 
velocity (in units of km s _1 ; not corrected for the heliocentric velocity). 




Fig. 5. — Intensity of the same lines depicted in Figure 3 as a function of the pixel position 
along a diagonal line that crosses region Rl, beginning at pixel (1,2), passing through region 
R3 (pixel 3,4) and finishing at pixel (11,12). Each type of line shows the intensity of a 
different spectral line: the solid line shows [N II] A5754; the short dashed line, the intensity 
of [N II] A6548; the dot-dashed line, Ha; the triple dotted-dashed line, [N II] A6583; the long 
dashed line, He I A7065 and, finally, the dotted line shows [Ar III] A7135. In order to show 
all of the intensities in the same plot, we have divided the Ha intensity by a factor of 10. 
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Fig. 6. — The flux-weighted rms width (top) and skewnes (bottom) moments of the intensity 
vs. velocity distribution, as given by equations (2) and (3), respectively, for the Ha (left) 
and [Ar III] A7135.8 (right) lines. The grey scale on the right side of each map is in unit of km 
s _1 . We note that both integrals show positive values and peak in a SE region, indicating 
the presence of a redshifted component towards the SE direction. The velocities are not 
corrected by the systemic radial velocity, and the integration limits in all maps are from 
< v ra d < 200 km s -1 . The white pixels have values close to zero. 
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Fig. 7. — a) Background and continuum subtracted data (full line) together with the Gaus- 
sian fit (crosses) for the Ha line at position (x,y)—(6,7) (see Figure 1). Top — left : Line 
profile and the three Gaussian fit. Top — right : Main, low velocity line profile and its Gaus- 
sian fit. Bottom — left : The data profile subtracted from the main fitted profile, showing 
evidence for a blue- and a red-shifted component. Bottom — right : The data profile sub- 
tracted from the three-component Gaussian fit. The residual profile is a random fluctuation 
around a flux of 0.5xl0 -15 erg cm~ 2 s _1 . These fluctuations corresponds to 10% of the flux 
for the blueshifted component, 0.5% of the flux for the main component and 10% of the flux 
for the redshifted component, b) The same but for the [Ar III] A7135 line. 
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Fig. 8. — The intensity of the main, low velocity peak (in units of 1 x 10 15 erg cm 2 
s _1 px" 1 ) as a function of position for the [N II]A5754 (top-left), [N II] A6548 (top-center), 
Ha (top-right), [N II]A6583 (bottom-left), Hel A6678 (bottom-center) and [Ar III]A7135 
(bottom-right) lines. Almost all of the lines peak at the proplyd position, and three of them, 
namely the [N II]A6548 (top-centre), Ha (top-right) and [N II] A6583 (bottom-left lines), have 
a secondary, less intense peak in region R3 (see Figure 1). 
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Fig. 9. — a) The intensity, in units of 1 x 10 erg cm s px , as a function of position 
in region Rl. b) The velocity, in km s _1 for the redshifted Gaussian fit to the Ha profiles, 
c) The jet velocity variation as a function of distance from the proplyd (see the text for a 
discussion). The velocities take into account the systemic radial velocity of v & ~ 26 km s -1 . 
The error associated with the radial velocity measurement is of the order of 2.5 km s _1 (see 
the text for a discussion). 
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Fig. 10.— a) N =([NII]A6548 + [NII]A6583)/[NII]A5754 ratio map of region Rl (see Figure 
1) of the observed field, b) A [T e ,n e ] diagnostic diagram showing the locci corresponding 
to different values of the nitrogen line ratio, c) The density for an assumed temperature of 
T = 10 4 K. d) The density for an assumed temperature of T = 1.5 x 10 K. 



